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We investigate the distribution of the supercurrent through a chaotic quan- 
tum dot which is strongly coupled to two superconductors when the Thouless 
energy is large compared to the superconducting energy gap. The distribu- 
tion function of the critical currents p{I c ) is known to be Gaussian in the 
limit of large channel number, N —* oo. For N = 1, we present an ana- 
lytical low-temperature expression for this distribution function, valid both 
in the presence and in the absence of time-reversal symmetry. It connects 
directly the distribution of transmission coefficients to the distribution of crit- 
ical currents. The case of arbitrary channel number (N > 2) is discussed 
numerically, and for small critical currents analytically. 

PACS numbers: 73.20.Dx, 74.50. +r, 74.60. Jg, 74.80.Fp, 73.20.Dx, 85.25. Cp. 



1. INTRODUCTION 

Electronic transport through mesoscopic structures with system size L 
smaller than the electron phase relaxation length L^, has become a unique 
"playground" to study quantum interference effects Jiro An intriguing phe- 
nomenon which convincingly presents the essence of quantum interference 
in diffusive electronic transport are the universal conductance fluctuations 
(UCF). Their size is of the order of e 2 /h, independent of sample geometry 
or the details of the realization of disorder, provided that the electronic sys- 
tem stays in the diffusive transport regime. UCF were first considered in a 
diagrammatic aDproach,Q'B and then soon discussed within random-matrix 
theory (RMT)!@ which ImryEI had proposed to apply for this class of prob- 
lems. UCF were observed as magneto-finger-prints in mesoscopic wires (see 
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the measurement of conductance fluctuations in gold wires by Washburn and 
Webb!). 

The superconducting analogue of UCF in normal conducting wires are 
supercurrent fluctuations in chaotic SNS-Josephson point contacts, as put 
forward by Beenakker. 00 The sample-to-sample fluctuations of the critical 
current I c depend only on the gap in the single particle spectrum: 

rms/ c ^ - <J C )2 = (1) 

with a numerical coefficient c of order unity. Here, A is the bulk gap, and 
the junction extension is assumed to be small compared to the coherence 
length. Again, the fluctuations are universal in the sense that they do not 
depend on junction parameters like, for example, the degree of disorder or 
the junction geometry. These universal supercurrent fluctuations still have 
to be verified experimentally. 

The realization depends on several conditions. First of all, the dynamics 
of the junction should be chaotic either due to impurities or boundary scat- 
tering. For example, for a diffusive junction, the junction length L (which 
is to be taken as the separation of SN interfaces) has to be much larger 
than the mean free path I in the normal-conducting part of the junction, 
I -C L. Secondly, to ensure universality, the superconducting energy has to 
be the only low energy scale. This amounts to the criterion that the Thou- 
less energy E c = h/r^ is much larger than the maximum gap in the energy 
spectrum, E c 3> A, where is the mean dwell time of electrons or holes in 
the junction. For a diffusive junction, this condition of a short dwell-time 
is fulfilled if L is much smaller than the coherence length, L <C £, where 
6 = (CoO 1 ^ 2 i s the diffusive coherence length, £ Q = hv F /(irA), and v F is the 
Fermi velocity. Moreover, the number N of transverse modes at the Fermi 
level which propagate through the junction, has to be large to allow an ex- 
pansion in the "dimensionless inverse conductance" (2e 2 /h)/(G) and thus 
to nrove that the numerical coefficient c in Eq. (|l|) is indeed independent of 

Physically, (a) multiple Andreev reflection at the interfaces)^ and (b) 
the statistics of the transmission eigenvalues lead to the universal fluctua- 
tions, Eq. (|]). Thereby, (a) is reponsible for the formation of current carrying 
bound (Andreev-) states in the junction. In the short dwell-time regime, only 
one Andreev state per channel exists. As a consequence, the supercurrent 
/(</>) is a linear statistic^ on the normal-region transmission coefficients T n 
(1 < n < N), independent of all junction parameters. However, the criti- 
cal current I c = m&xl(<f)) = I((j) c ) is in general not a linear statistic on T n 
since the phase cj> c at which the supercurrent is maximal depends on all the 
transmission coefficients. Beenakker showed that I(4> c ) approximately obeys 



RMT: Distribution of Mesoscopic Supercurrents 



a linear statistic to within a correction of order 1 /N, where N is the number 
of channels. Prom the statistical properties of the transmission eigenvalues 
T n it is known that the fluctuations of a linear statistic are of order unity 
and thus are universal, Eq. ([!]). 

In this work we consider a Josephson junction where a chaotic quan- 
tum dot is strongly coupled to two superconductors. It is known that the 
scattering matrix of the normal conducting dot in this case belongs to the 
"circular unitary ensemble" (CUE) of RMT.Ej The case of a weakly cou- 
pled junction was considered elsewhereJll We are especially interested in the 
fluctuations of the supercurrent if the number of conducting channels in the 
normal chaotic region, N, is small. This implies that the critical current 
does not obey a linear statistic. Nevertheless, we expect that the sample-to- 
sample fluctuations are still determined by the relation Eq. (|l]). However, 
the prefactor c will now be weakly dependent on the number of channels N. 
More relevant is the question of how the full distribution of critical currents, 
p{I c ), will behave when we transit from a situation with no level repulsion 
[N = 1) to the scenarios with level repulsion [N > 2). Certainly, the two 
situations differ qualitatively. For N —* oo the distribution is a Gaussian.0 
For a single-channel junction, the literature does not provide a numerical or 
analytical result, to our knowledge. The missing level repulsion will move 
the center of the distribution to a lower value of I c , the details of which 
should depend on the chosen ensemble. Finally, the case of a small num- 
ber of channels is intriguing because it will interpolate between these two 
qualitatively incompatible situations — no level repulsion versus repulsion 
between many levels. 

Although a Josephson junction with these properties has not yet been 
realized, neither for a small nor for a large number of channels, the mi- 
crofabrication of point contacts with a small number of channels is already 
feasible (for example, see Refs. 15.19.2C); and findings like the conductance 
quantization, related to the Friedel sum rulej^il open new possibilities for 
the research of correlation effects in nano-scale structures. Chaotic point 
contacts are now within reach of experiments. The investigation of the in- 
terference of diffusion and Coulomb interaction should soon be possible, and 
we expect new exciting physics to emerge in these systems. 



2. MODEL 

In the recent 15 years, RMT has become a tool with which electronic 
transport through disordered mesoscopic systems can be studied very effi- 
ciently. It builds on the statistical independence of the elements of, e.g. the 
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scattering matrix and directly implements the symmetry of this matrix. As 
an excellent review by BeenakkertS on this topic is available, we will only 
motivate some of the basic results in this section. We are going to apply 
RMT to the normal region of the SNS contact. The model which we use 
for the SNS contact was suggested in Ref. 11. However, we will extend the 
previous considerations to heterogeneous contacts. 



2.1. Random Matrix Theory for the Chaotic Region 

The critical current through an SNS-junction is determined by a specific 
realization of a random potential in the normal region which is the origin of 
the mesoscopic fluctuations. Interaction effects between the electrons in the 
normal metal are not taken into account, neither is spin-orbit coupling. 

Concerning the interface between the superconducting and the normal 
regions, two scenarios have to be distinguished: the case of weak coupling, 
realized for instance by tunneling barrierst^ and the case of strong coupling 
with ideal, impurity-free leads at z\ and (see Fig. [l]). In the latter case, the 
scattering matrix S is a unitary random matrix with eigenvalues uniformly 
distributed on the complex unit circle,tl'ii and the corresponding ensemble is 
referred to as CUE. In this article we will consider only the circular ensemble, 
that is, we stay in the strong coupling regime. Furthermore, we assume the 
Thouless energy to be much larger than the gap, E c 3> A. In this limit, which 



was named "short dwell-time regime" in Ref. 16, the energy dependence of 
the scattering matrix can be neglected. 

On the other hand, the dwell-time Td of electrons and holes in the chaotic 
region — with energies close to the Fermi level - has to be sufficiently large 
so that RMT is applicable. The dwell-time Td should be larger than the time 
scale T erg at which ergodicity sets in, where h/r eT g = hv F mm{l, L}/L 2 . 

The scattering matrix S has the block structure 

(2) 



r if 
t r' 



with the reflection matrices r, r' and the transmission matrices t, t', all of 
which are N x N matrices. The number of propagating modes in each of 
the two leads is N. Current conservation implies that the scattering matrix 
is unitary: = S* , and a "polar decomposition" is possible!? 

s = ( o v){ ^VT^ J^T ) ( V' ) ' (3) 



where a N x N diagonal matrix T = di&g(Ti,22, . . . , Tjv) has been intro- 
duced. The transmission eigenvalues T n are the eigenvalues of the matrix 
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superconductor 1 normal conducting superconductor 2 




Fig. 1. Junction geometry and definition of scattering amplitudes for parti- 
cles (e) and holes (h). 

product ttf = t't^ . The four unitary N x N matrices U,U' and V, V are 
to be interpreted as "phase factors". In the case of time-reversal symmetry 
the scattering matrix is unitary and symmetric, S = S T , which implies that 
also V = U T and V = V T hold ("circular orthogonal ensemble" (COE)). 

For statistical averaging we need the Haar measure fi which can be 
determined for the above polar decomposition asllil 

dfi(S) = J({T n }) J] dn{U a ) I] dTi . (4) 

a i 

The Jacobian is the probability distribution P({T n }) of the transmission 
eigenvalues (up to a normalization constant cn); it is given by 

J({T n }) = II \Tn-T m fl[T k 1+m = c N P({T n }) (5) 

n<m k 

Here, f3 is the ensemble parameter — which numbers the degrees of freedom 
of the matrix elements of S: 

(5 = 1 COE system with time-reversal symmetry 

(3 = 2 CUE no time-reversal symmetry (6) 

for the two cases considered in this paper. The set of independent unitary 
matrices in the Haar measure (||) is {U a } = {U, V,U',V'} for (5 = 2, and 
{U a } = {U,V} for (3 = 1. For the purpose of interpreting the probability 
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distribution (||) it has become popular to write it in the form of a Gibbs 
distribution 



P({T n }) = — exp ((3 ln \ T n ~ T m \ -(3^y(Tk)) 

CN \ n<m k J 



(7) 



with the "one-particle" potential V(Tk) = (1/(3 — 1/2) In and a "two- 
particle" logarithmic repulsion. 



2.2. Scattering Matrix of the SNS Junction 

Here we follow the approach pioneered by Beenakker,0 and extend it in 
order to include the case of heterogeneous junctions.0 The details of how to 
set up the Bogoliubov-de Gennes (BdG) equation for this type of junction 
are found in the cited paper and in Ref. [25|, and will not be repeated here. 

The S-matrix of the complete junction, S SNS , connects in- with out- 
going amplitudes of the electronic wave function in the superconductors. 
The components of the wave function are defined in Fig. [I]; for example, we 
have for the outgoing amplitudes N components a~ for electron-like quasi- 
particles running from the interface z\ to the left and N components d+ 
for electron- like quasiparticles running from the interface Z2 to the right. 
Andreev reflection mixes particle with hole components, and therefore the 
corresponding hole-like quasiparticles with N components and have 
to be included in the scattering process. The AN x 4A r -matrix S SNS is now 
defined through 

/ a7 \ ( at \ 



4 
V 4 



V 4 J 



(8) 



The matrix S SNS is actually constructed from more basic 4A^x4A^ S-matrices, 
S N and S A . The matrix S N is formed by the scattering matrices for electrons 
and holes in the normal, disordered region: 



S N (£) 



S(E) 
S*(-E) 



(9) 



where E is the energy relative to the Fermi energy. The matrix S A , defined 
by 



e -i((p(£0+#) 



Sa(E) 




(10) 
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is the scattering matrix for Andreev reflections at both NS interfaces. The 
off-diagonal structure originates in the fact that Andreev reflection trans- 
forms an electron into a hole. The transformation is supplemented by two 
phase shifts: one from the penetration of the wave function into the super- 
conductors 

\ U ^2 ljV / 

where 1/v is the N x N unit matrix, and 

/ 

(fj = —i In 



E 




IT ~ 1 ■ (12) 



The Aj are the respective Bogoliubov quasiparticle energy gaps of the two 
superconductors. We use rigid boundary conditions with step-like gaps at 
the SN-interfaces. This is approximately valid for the considered case that 
the extension of the junction is much smaller than the coherence length in the 
bulk superconductor.!! The second contribution to the phase shift in (10) is 



the phase 4>i^ acquired from the pair potential in the superconductors which 
has the matrix form 

*=f 01 o" A.V (13) 



>2 J- TV 



With these basic elements the scattering matrix for multiple Andreev re- 
flection at the interfaces, in combination with the intermediate propagation 
through the disordered region, takes the form0 

S n (E)[l m -Sl(E)S A (E)} 

( -i( v -*)/2 o \ 

[l^-S^S^)]" 1 Q e -^+*)/ 2 (14) 




where I47V is the 4N x 4N unit matrix. Obviously, the ipj defined in (12) 
are real for E < min{Ai, A2} = A min . In this case, the matrix S A is unitary: 
Andreev reflection within the region between the two interfaces is current 
conserving. For \E\ > max{Ai,A2} = A max , the matrix S A is hermitian, 
and only then is S SNS a unitary matrix. 



2.3. Spectrum of States of a SNS-Junction 



Corresponding to the range of energies at which we consider the electron 
propagation through the junction, we encounter three different regimes: 
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1. \E\ > A max , in which case we consider a contiuum of particle and hole 
scattering states in each superconductor, 

2. A max > E > A min , which is the intermediate regime with quasiparticle 
states extending over the superconductor with the lower energy gap 
but they decay exponentially in the superconductor with the larger 
gap, 

3. \E\ < A min , in which case the BdG wave functions decay exponentially 
in both superconductors. Bound states, Andreev states, which are 
localized in the normal region, generate a discrete spectrum. For a 
finite phase difference, 4> = 4>2— $1, Andreev states support an electrical 
current. Thereby, the coherent particle- hole excitations in the junction 
are converted into a supercurrent in the bulk superconductors. 



2.3.1. Continuum states 

In order to obtain the spectrum in this regime, one makes use of the 
"Wigner time-delay matrix'Sil S^^S, the trace of which is the density of 
states, p(E): 

1 f . d 1 Id 

p(E) = tr SO = ^hidetS. (15) 

FV 1 2m \ dE J 2m dE y ' 

As the scattering matrix of the junction S SNS is unitary, one may evaluate the 
relation ( |T5| ) with S = S SNS and find, respecting the hermiticity of S A (E), 

P(E) = ~^lnZ(E)+p (E) (16) 
with Z(E) a complex function defined through 

Z(E) = det(l 4JV - S A (E)S N (E)) . (17) 

The additive term 

» (E) = i4 lndets » < 18) 



in the density of states fllq ) is independent of the phase difference (f> of the 
two superconductors and will not contribute to the Josephson current. 



2.3.2. Intermediate regime 



The determination of the density of states in this energy regime does 
not proceed parallel to that of p(E) with extended quasiparticle states in 
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both superconductors. The scattering matrix of the junction is not unitary 
for A max > E > A min since the quasiparticle states decay exponentially in 
the superconductor with the larger energy gap — we assume this to be the 
left superconductor in Fig. |l[ However, we can define a reduced scattering 
matrix S int which connects only the incoming amplitudes from the right 
superconductor to the outgoing amplitudes of the right superconductor 



(19) 



and the dimension of this scattering matrix is half that of S S ns- Since S lnt is 
unitary we may use relation (|l^) for S = S int . After some tedious algebra 
we arrive again at 

p{E) = -±^lnZ(E)+p (E) (20) 
with a ^-independent additive term p~Q 

As already indicated, the scattering matrix for Andreev scattering S A (E) in 
Z{E) is hermitian for the continuum states but not so for the intermediate 



states, and expressions (16) and (EG) are indeed quite distinct. 



2.3.3. Bound states 



Since the Andreev states are evanescent in the superconductors, the 
matrix S A (-E') is now unitary. It obeys the simple relation for bound states 



K 
V c t ) 



S A (E)S N (E) 



K 
\ 4 J 



(22) 



Solutions exist if Z(E) = is fulfilled. In this regime the density of states 
can therefore be written as 



P(E) 



dZ 



i 



dE 



l^-Q]nZ(E + iO) 

7T dE 



it dE Z(E + iO) 



(23) 
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2.4. Josephson Current 

The stationary current I through a SNS-j unction is driven by the phase 
difference of the superconductors. It can be obtained from a free energy F 
through the well-known thermodynamic relation 

j-^ (24 ) 



As elaborated in the literature (for example, see Ref. (|Tl|)) the current is 
given in terms of the density of quasiparticle states 

oo 

I = 2k B r ( dE In [2 cosh (E/2k B r)] —p{E) . (25) 
Ji J deb 

o 

After an integration by parts the current can be expressed as 

oo 

I = -| f dEtanh[E/2k B T]j(E + i0,4>) , (26) 
Tin J 
o 

where j is a kind of spectral current 

j(E,</>) = -^lnZ(E) 

= --^lndet(l 4 7v-S A (£)S N (£)). (27) 

dep 

This compact expression for the supercurrent is the basis for the calculation 
of the distribution of supercurrents through chaotic junctions in the following 
sections. 



3. SINGLE CHANNEL CASE 

Quantum transport of coherent particle-hole excitations through a sin- 
gle channel is rather particular. Level repulsion is missing by construc- 
tion. The transmission eigenvalues T are distributed according to a "single- 
particle" potential which is either logarithmically diverging for small T 
(P = 1) or constant ((3 = 2), P(T) = |t~ 1+/3 / 2 . This potential will control 
the distribution of the critical currents. It should be realized that this sce- 
nario applies also for junction geometries where (many) different channels 
are uncoupled. 

It is quite instructive to derive step by step the spectral current for 
this single channel case since the procedure conveys an understanding of the 
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validity and consequences of the approximations made. Here we present the 
result for the case of a scattering matrix which is energy independent, that is, 
the Thouless energy is much larger than A max , and there are no resonances 
at or in the gap. The following relation holds 

j(E,(ft) = -— 3?ln[cos(^i(£) + ip 2 (E))- Rcos^E)- <p 2 (E))- Teas fa], 

(28) 

which is derived in Appendix A where we also discuss an expansion in E/E F 
of the phases of the scattering matrix. In the above expression, we introduced 
R = 1 — T and </> e g- = (ft + A(ft, where A(ft is a phase which is generated from 
the phase factors U,U' and V,V in the polar decompostion (||). Since A<ft 
builds only on a difference of the respective phases, it is zero for the COE 
where U' = U T and V = V T hold, whereas it runs from to 2ir for the 
CUE. 

We now restrict our discussion to symmetric contacts with Aj = A2 = 
A. The energy eigenvalues for the Andreev states are straightforwardly 

found from the zeros of the logarithm in (^), E((ft) = ±Ayl — Tsin 2 ^. 



2 • 

The supercurrent (^6|)_is carried only by the Andreev states (the continuum 
contributions cancel) 



„x eA T sin (ft , . A / m 6 S , . 

2h ^l-Tsin 2 ! 2k B T \ 2 

Obviously, this expression recovers the well known sine-phi dependence for 
small transmission coefficients, I ((f)) = ^Gsin (ft, with the conductance G = 
2e 2 T/h. 

Finally, the current for a particular junction (p9|), characterized by a 
specific T, is averaged with the distribution function (0) for a single channel 
and we find for r = and (3 = 1 (see Fig. |2[ lowest curve) 

(^(</>)}coe = sin (ft ^ . Sin / for (ft G [-vr, vr] (30) 
n 8 sm d § 



For (3 = 2, the phase difference (ft in fl29|) has to be replaced by (ft e g. Inte- 
gration over dfj,(U a ) (see Eq. (Q)) corresponds to averaging over A(ft which 
results in a zero ensemble-averaged supercurrent for the CUE. 

Moreover, the variance of the supercurrent can be calculated analyti- 
cally. However we want to focus on a different issue: the critical current, that 
is, the maximum supercurrent is usually more accessible in experiments. In 
the single-channel case the structure of the supercurrent-phase relation ( p9[ ) 
is actually so elementary for zero temperature that it is possible to determine 
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CUE 


COE 


P(T) 


1 


1 




2(l-i c ) 


(1 " ic) 
^Ji c (2 - i c ) 




1 

3 


4 


var i c 


1 

18 


2 7T 2 

3 16 



Table 1. Disbribution of the critical currents for the CUE und COE; the 
average critical current is (i c ) = J 1 i c/ o(i c ), and the variance is defined as 
vari c = (i c 2 ) - (i c ) 2 . 

directly the distribution function for the critical current. The distribution 
function p(i c ) for the dimensionless critical current, i c = J c /^-, characterizes 
the ensemble in a unique way, as we will present below. A suitable definition 
of the distribution function for the critical current is 

i 

p(i c ) = J dT P(T) 5(i(T, MT)) ~ ic) • (31) 
o 

The critical current for a specific transmission coefficient, i(T,<p c ), is found 
from ( pgp (at r = 0) together with the condition di((p c ,T)/d4> c = 0. It 
yields an implicit equation for <f> c which reads T = T C (0 C ) = — cos <p c / sin 4 
This latter relation allows to express the critical current in terms of 4> c : it 
reads i(T c ((j> c ), C ) = 1 — (tan4r) -2 . Since the function T c ((j) c ) is monotonic 
and maps the interval [tt/2,tt] onto the interval [0,1] we may rewrite the 
integral (^) in terms of an integration over <\> c . Finally, making use of 
^c(0c) = ^c(2 — i c )i we find the fundamental zero temperature relation for 
the distribution of the critical current: 

p(i c ) = 2(1 -i c )P(i c (2 -i c )). (32) 

Again, i c is normalized to the maximum physical critical current eA/K, so 
it covers the interval [0, 1] . We would like to emphasize that this formula 
is valid in the presence as well as in the absence of time-reversal symmetry. 
For the single-channel case breaking of time-reversal symmetry results in a 
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shift of the I((f>) curve by a constant phase A^ which does not affect the 
critical current. 

The immediate consequences for the COE and CUE are summarized 
in Table ||. Relation (|32] ) allows to determine the type of ensemble if the 
distribution of critical currents is measured. The distribution function for the 
COE diverges as l/\/^c for i c — ► (which is also displayed in Fig. ||, N = 1), 
whose origin is the square root divergence of the distribution function P(T). 



4. MULTIPLE CHANNEL JUNCTIONS 

Only in the presence of time-reversal symmetry, (3 = 1, the expression 
(|28| ) for the spectral current can be generalized to multiple channels, 

N , 

j(E,(j>) = -Y^ ^ln[cos{^ 1 {E)+ip 2 (E))-R n cos(ip 1 (E)- t f2(E))-T n cos(t ) }. 

n=l " 

(33) 

In Appendix B the averaged critical supercurrent for the asymmetric Joseph- 
son junction is computed from this expression. For the symmetric junction, 
Ai = A2 = A, again only the Andreev-states contribute and the current is 
given as a linear statistic as obtained by Beenakker,lill 

K*. IH}) -ft tanh(^^/l - W f ) ■ (34) 

The average supercurrent can be determinedanalytically in the limit 
of large channel number N and zero temperatureS It is obtained with the 
density of transmission coefficients, p(T) = N/ir [T(l — T)] _1 , and is non- 
sinosoidal due to its dependence on the hypergeometric function 2-F1, 

W))co B = N ^ 2^1(^,2, sin 2 |) sin <j> for N 00. (35) 

This expression, valid for large N, is to be compared to results for 
finite N. For this purpose, the average supercurrent has been evaluated 
numerically for N = 1 to N = 5. Although the curves do not scale with 
N — the maxima are slightly shifted with increasing N towards higher values 
of <fi (see Fig. |2|a) — the overall shape of the supercurrent-phase relation is 
remarkably similar for N = 1 to N = 00. The same holds, to a lesser 
extent, for the root-mean-square (rms) fluctuations (Fig. |||b). For N = 00, 
we included the result obtained by Beenakker,E2l which is slowly approached 
by the finite- N results. He argued that corrections in 1/N should move the 
finite value of rms I^ =00 {(p = tt) = eA/irh to zero. 
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Fig. 2. Average supercurrent for the COE, scaled with NeA/h, as a function 
of the phase difference <ft across the junction (a), and the rms fluctuations of 
the supercurrent, also for the COE (b). 



Much more pronounced is the channel-number dependence of the distri- 
bution function for the critical current. The divergence of the distribution 
function in the COE is lifted in the multiple-channel situation due to the 
interaction between the transmission eigenvalues, that is, due to level repul- 
sion. It is not possible to derive a closed expression for the distribution in 
the case of two or more channels. However, an analytical expression can be 
found which controls the limit i c — > 0. The proof exploits again the simplic- 
ity of the current-phase relation (^4|) where we restrict the evaluation to zero 
temperature. The limit i c — » is approached when all T n — > 0. In this case, 
we neglect the square-root denominator, which results in <p c — > tt/2. As ex- 
pected, the critical current is in this limit proportional to the conductance g, 
i.e. the sum over all transmission eigenvalues, i(<j> c , {Ti}) = \ J2n=i T n = \g- 
Therefore, for small i c , the distribution of the critical current behaves in the 
same manner as the distribution of the conductance does.il 

Similar to the single-channel situation, the distribution function for the 
critical current is defined as 

P(ic) = c^ 1 J] [ dT l II \Tn-T m \Y[T~ 1/2 5(i(M{Ti}),{Ti})-ic) 

l Jo n<m k 

(36) 

After rescaling all transmission coefficients in ( |36|) with i c one obtains the 
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I c [NeA/h] 

Fig. 3. Distribution function for the critical current in the COE for various 
channel numbers N. The maximum value reached by the averaged {I ((ft)) 
(Fig. ||a) curve has been marked in the graph with the corresponding channel 
number by a black dot. The inset shows the distribution for N = 2 and its 
derivative. 

limiting small-current behavior 

N 2 

p(i c -> 0) ~ i c a with a = — - 1 (37) 

For a single channel, N = 1, this agrees with the result of Table [l]. In 
particular for two channels one obtains 

p(i c 0) ~ i c for N = 2. (38) 

Level repulsion is responsible for the vanishing distribution function in the 
limit of zero critical current. 

The distribution functions for the full range of i c are generated numeri- 
cally via a Monte Carlo integration. The results are displayed in Fig. [|. The 
numerical error is of the order of the line width. The distribution for N = 2 
again is shown in the inset together with its derivative. The derivative indi- 
cates clearly the feature at I c = eA/h. For larger channel number N, the dis- 
tribution functions approach a Gaussian distribution. In the limit N — > do, 
the Gaussian is characterized by (J c ) = 0.35 NeA/h and rms/ c = 0.30 eA/h, 
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and the critical phase, at which the supercurrent attains its maximum value, 
is 4> c = 2.04. For N = 10 we are already close to this limit, with first mo- 
ments (J c ) = 0.33 ± 0.01 NeA/h and rms/ c = 0.29 ± 0.05 eA/ft. 



5. SUMMARY 

The distribution function for the dimensionless critical current i c = 
I c /(eA/h) throirgh a chaotic quantum dot is, for a large number of channels 
N, a GaussianJUl Its position is the average critical current of order TV and 
its width are the root-mean-square fluctuations of i c of order unity. With 
decreasing channel number these two first moments of the distribution func- 
tion are still of order N and 1, respectively, however the prefactors change. 
Both, the average value of the critical current, normalized to N, and the rms 
fluctuations of the critical current decrease. This is not unexpected since the 
"level repulsion" between the transmission eigenvalues of different channels 
declines. 

In the extreme case of a purely one-dimensional normal conducting con- 
tact, that is a single-channel contact, level repulsion is entirely absent and 
for the COE the distribution function displays a square-root divergence. 
This originates in the singular behavior of the distribution function for the 
transmission eigenvalues for the COE. Moreover, a general relation was es- 
tablished which connects the distribution of transmission eigenvalues to the 
distribution function for the critical supercurrent in the single channel situ- 
ation. For multiple channels we provided a small current expansion of this 
distribution function which relates the number of channels to the exponent 
of the algebraic small-current behavior. With these relations it is in princi- 
ple possible to identify the ensemble by a measurement of the distribution 
of critical currents. 
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APPENDIX A 

Here, we elaborate on the determination of the spectral current j(E, <f>), 
Eq. (^), for the single channel situation. An understanding of the energy 
dependence in this most simple case is useful for further discussions. The 
typical energy scale for variations of the scattering matrix elements is the 
Thouless energy. 

The scattering matrix S for the normal region of the contact is a 2 x 2 
matrix. In the polar decomposition (|3|) it is characterized by four phases 
Ol, Or, XLi XR an d the transmission coefficient T of this single channel: 

„ = ( e idL \ / -VT^T VT 

In general, these five parameters of the scattering matrix are energy depen- 
dent. We define the new phase parameters = Ol + Or, AO = 6l — Or, 
X = XL + XR an d Ax = XL — XR an d find for the spectral current 

j(E,<p) = -^ln[cos(ip 1 (E)+ME)-l(O(E)-0(-E)+x(E)-x(-E))) 

-R(E) cos(^i(£) - <f 2 (E) - ^(0(E) - 0(-E) - X (E) + x(-E))) 

-f(E) cos(0 + \{A0{E) + AO(-E) - A X (E) - A x (-E)))} (40) 

with the particle-hole symmetric values for transmission and reflection coef- 
ficients 

f(E) = ^T{E)T{-E) 

R(E) = y/(l-T(E))(l-T(-E)). (41) 

If the matrix S is only weakly energy dependent (in particular no reso- 
nances should exist in an energy range of 2 A around the Fermi energy Ep) 
the above expression can be expanded: according to the Andreev approx- 
imation the phases are taken to first order in E/Ep and the transmission 
coefficient in zeroth order. This results in the low-energy expression for the 
spectral current 

j(E,<t>) = -±Qln[cos( l p l (E) + ^2(E)-(0'(0) + X '(0))E) 
acp 

-Rcosi^E) - ME) ~ (0'(0) - X '(0))E) 

-Tcos(0 + A0(O) - Ax(0))] (42) 

where the prime on the phases denotes a derivative with respect to energy, 
and R = 1 - T. 
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Fig. 4. Critical current for large N of an asymmetric Josephson contact 
averaged over the COE and of a ballistic contact for different ratios 5 = 

^min/^max ■ 



Only in the presence of time-reversal symmetry the phases A6(0) and 
Ax(0) vanish. However, in the absence of time-reversal symmetry the whole 
I ((f)) curve is shifted by A0(O) — A%(0). The mean supercurrent at a fixed 
phase (j) then requires averaging over the constant phase shift which results 
in a vanishing supercurrent. On the other hand, a measurement of only the 
critical current is not affected by this phase shift and therefore expression 
(|32| ) is valid irrespective of the presence or absence of time-reversal symme- 
try. 

Assuming the transmission coefficients to be distributed according to 
the COE, the averaged spectral current fll^ ) reads in the presence of time- 
reversal symmetry 



where 



COE — 



8111 ( 



[ cos (f) - cos(9?i - tf2 — (0' - x')E) 

C0S(<£1 — If 2 — - X')E) - COS(j) 



(1 



arctan y/a. 



cos(v9i +ip 2 - (0' + x')E) - cos(^i -ip 2 - {0' - x')E) 



The derivatives of the phases are taken at E = 0. The averaged spectral 
current (fEj) can be considered as a basis for a more detailed evaluation 
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respecting a weak energy dependence of the scattering matrix. 

If we neglect the energy dependence of the scattering matrix altogether, 
we arrive at Eq. (28), which was the basic equation for the evaluation in the 
section on single channel junctions. 



APPENDIX B 



In the limit of a large number of channels N the expression fl33|) for the 
spectral current of an asymmetric Josephson junction can be averaged over 
the COE with the distribution p(T) = N/tt [T(l - T)]" 1 . We get 



</>)) = 



iVsin< 



cos (j) — cos (ipi — tp2 ) 



I COs((^i — if2) — COS (j) 

cos(ipi+if 2 ) - cos(</?i— <p 2 ) 




The averaged critical current (I c ) for large N can now be evaluated numeri- 
cally and is shown in Fig. Q versus the asymmetry parameter 5 = A min /A max . 
For 5 = 1 it starts at the value (J c ) = 0.35 NeA/h for the symmetric junc- 
tion and decreases with decreasing 5. For comparison the values of I c for 
the ballistic junction with N open channels are also shown. 



REFERENCES 

1. B.L. Al'tshuler, P.A. Lee, and R.A. Webb, eds., Mesoscopic Phenomena in 
Solids (North-Holland, Amsterdam, 1991), and references therein. 

2. Y. Imry, Introduction to Mesoscopic Physics, (Oxford University Press, 1997). 

3. For an extended discussion of mesoscopic Josephson effects see: Proceedings 
of the NATO Advanced Research Workshop on Mesoscopic Superconductiv- 
ity, eds. F.W.J. Hekking, G. Schon and D.V. Averin (1994, North-Holland, 
Amsterdam, or Physica 203B, Nos. 3 & 4). 

4. B.L. Al'tshuler, Pis'ma Zh. Eksp. Teor. Fiz. 41, 530 (1985) [JETP Lett. 41, 
648 (1985)]. 

5. P.A. Lee and A.D. Stone, Phys. Rev. Lett. 55, 1622 (1985). 

6. Muttalib, Pichard, and Stone, Phys. Rev. Lett. 59, 2475 (1987). 

7. P.A. Mello, P. Pereyra, and N. Kumar, Ann. Phys. (N.Y.) 181, 290 (1988). 

8. Y. Imry, Europhys. Lett. 1, 249 (1986). 

9. S. Washburn and R.A. Webb, Adv. Phys. 35, 375 (1986). 

10. C.W.J. Beenakker, Phys. Rev. Lett. 67, 3836 (1991) and Phys. Rev. Lett. 68, 
1442 (E) (1992). 

11. C.W.J. Beenakker, Three "universal" mesoscopic Josephson effects in: 
Proc. 14th Taniguchi Int. Symp. on 'Physics of Mesoscopic Systems', eds. 
H. Fukuyama and T. Ando (Springer, Berlin, 1992). 



M. Garst and T. Kopp 



12. A heterogeneous SNS-junction consists of two superconductors with different 
values of the superconducting gap, Aj 7^ A 2 , whereas for a homogeneous 
junction Ai = A 2 . 

13. A. F. Andreev, Zh. Eksp. Teor. Fiz. 46, 1823 (1964); 49, 655 (1965) [Sov. 
Phys. JETP 19, 1228 (1964); 22, 455 (1966)]. 

14. A quantity A — ^ n a(T n ) is denoted a linear statistic on the transmission 
coefficients because products of T n are absent. 

15. C.W.J. Beenakker, Rev. Mod. Phys. 69, 3 (1997). 

16. P.W. Brouwer and C.W.J. Beenakker, Chaos, Solitons & Fractals 8, 1249 
(1997). 

17. H.D. Politzer, Phys. Rev. B 40, 11917 (1989). 

18. E. Scheer, W. Belzig, Y. Naveh, M.H. Devoret, D. Esteve, and C. Urbina, 
Phys. Rev. Lett. 86, 284 (2001). 

19. M.F. Goffman, R. Cron, A. Levy Yeyati, M.H. Devoret, D. Esteve, and 
C. Urbina, Phys. Rev. Lett. 85, 170 (2000). 

20. M.C. Koops, G.V. van Duyneveldt and R. de Bruyn Ouboter, Phys. Rev. Lett. 
77, 2542 (1996); mechanically controllable break junctions, which are indeed 
atomic-size quantum point contacts, were developed first by C.J. Mullcr, J.M. 
van Ruitenbeek and L.J. de Jongh, Physica 191C, 485 (1992). 

21. S. Kirchner, J. Kroha, and E. Scheer, Generalized conductance sum rule in 
atomic break junctions in: Proc. of the NATO Advanced Research Workshop 
'Size dependent magnetic scattering', Pecs, H ungary, May 28 - J une 1, 2000, 
in press (Kluwer Academic Publishers, 2001); 3ond-mat /00 1 1 p"^ . 

22. Th. Martin and R. Landauer, Phys. Rev. B 45, 1742 (1992). 

23. H.U. Baranger and P.A. Mello, Phys. Rev. Lett. 73, 142 (1994). 

24. R.A. Jalabcrt, J.-L. Pichard, and C.W.J. Beenakker, Europhys. Lett. 27, 255 
(1994). 

25. K. Bottcher and T. Kopp, Phys. Rev. B 55, 11670 (1997). 

26. K.K. Likharev, Rev. Mod. Phys. 51, 101 (1979). 

27. M.L. Goldberger and KM. Watson, Collision Theory (John Wiley & Sons, 
1964). 

28. E. Doron and U. Smilansky, Phys. Rev. Lett. 68, 1255 (1992). 

29. W. Haberkorn, H. Knauer, and J. Richter, Phys. Status Solidi A 47, K161 
(1978). 

30. C.W.J. Beenakker, Phys. Rev. B 47, 15763 (1993). 



